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We examine a dark energy model where a scalar unparticle degree of freedom plays the role of 
quintessence. In particular, we study a model where the unparticle degree of freedom has a standard 
kinetic term and a simple mass potential, the evolution is slowly rolling and the field value is of 
the order of the unparticle energy scale (X u ). We study how the evolution of w depends on the 
parameters B (a function of the unparticle scaling dimension d u ), the initial value of the field fa 
(or equivalently, \ u ) and the present matter density Q m o- We use observational data from Type 
la supernovae, baryon acoustic oscillations and the cosmic microwave background to constrain the 
model parameters and find that these models are not ruled out by the observational data. From a 
theoretical point of view, an unparticle dark energy model is very attractive, since unparticles (being 
bound states of fundamental fermions) are protected from radiative corrections. Further, coupling 
of unparticles to the standard model fields can be arbitrarily suppressed by raising the fundamental 
energy scale Mf , making the unparticle dark energy model free of most of the problems that plague 
conventional scalar field quintessence models. 



I. INTRODUCTION 

Cosmological data from a wide range of sources includ- 
ing type la supernovae [H, 0, HI , the cosmic microwave 
background [J, baryon acoustic oscillation s [a, Bl , cluster 
gas fractions [3, Q and gamma ray bursts [9l.ll0| seem to 
indicate that at least 70% of the energy density in the uni- 
verse is in the form of an exotic, negative-pressure com- 
ponent, called dark energy. (See Ref. [ld| for a recent re- 
view) . It is usually assumed that a scalar field with a pos- 
tulated potential may yield the equation of state needed 
to drive the acceleration. While it is easy to come up with 
such a potential (see e.g. [13, El, El, El, EH ) , the under- 
lying problem is that we have not discovered a single fun- 
damental scalar field so far. Fundamental scalar fields, 
being unprotected from radiative corrections, are very 
problematic since they require a significant fine tuning in 
theory (see e.g. [13, EH)- Couplings of the quintessence 
field with other standard mod el p articles, in particular 
fermions, must be suppressed [lj| HI, HH, HH (however, 
see [H, [3, [2f| for a discussion of PNGB quintessence 
models which are radiatively stable and see [H, HH for 
a quintessence model which can couple strongly enough 
with standard model particles to be detected in colliders). 
It is therefore instructive to search for other degrees of 
freedom which may effectively play a role similar to the 
scalar field. 

Recently, the existence of a new scale invariant sector 
very weakly coupled to the standard model was postu- 
lated [13, El, HI]. The fundamental energy scale, Mp, of 
this sector is perhaps far beyond the reach of today's or 
near-future accelerators. However, the existence of such 
a sector may affect low energy phenomenology. The effec- 
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tive low energy field theory which describes these effects 
is often called unparticle physics since these new degrees 
of freedom would not behave as ordinary particles. For 
example, their scaling dimension does not have to be an 
integer or half an integer. Because of the scale invariancc, 
the fundamental particles arc masslcss. An unparticle is 
a composite state of the fundamental massless particles, 
and couples to the standard model particles through a 
heavy mediator. This heavy mediator has a mass of the 
order of Mp. Thus, interactions between the unparti- 
cles and the standard model particles is suppressed by 
powers of Mp. Another important characteristic of un- 
particles is that their mass can take continuous values 
[301 ] - For example, the mass of an unparticle with en- 
ergy uj can take all the values from zero to w. A lot 
of work [3ll [H, \M, [H, [H recently focused on the 
new collider signals for unparticle physics. The unpar- 
ticles could also play the role in gravity and cosmology 

[33, m, m, m m m, m, m, m, m, S3, m • Black hole 

Hawking radiation of unparticles was studied in [4{| . 

In this paper we will try to answer the question 
whether scalar unparticle degrees of freedom can play the 
role of dark energy. An unparticle dark energy model is 
interesting for a number of reasons. In the simple model 
of unparticles constructed in [28[, the fundamental de- 
grees of freedom at high energies are fermions. Unparti- 
cles are bound states of these fermions, but at energies 
below a certain "phase transition scale" A„ they appear 
as the fundamental degrees of freedom. This situation is 
analogous to quarks and mesons in QCD. While funda- 
mental scalar fields are not protected from radiative cor- 
rections, unparticles are, since (like in technicolor mod- 
els) the fundamental fields are fermions whose mass is 
protected. Also, coupling of unparticles to the standard 
model fields can be arbitrarily suppressed by raising the 
fundamental energy scale Mp. Therefore, the unparticle 
dark energy model is free of most of the problems that 
plague fundamental scalar field quintessence models. 

In exploring this idea further, we will study in detail 
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an unparticle dark energy model with </> 2 quintessence 
potential. Since the effective scalar unparticle degree of 
freedom <j> is quite different from the fundamental scalar 
(i.e. different scaling dimension, continuous mass etc.) 
the phenomenology may be quite different from the stan- 
dard one. We will derive the effective functional depen- 
dence of the energy density of the universe on </> and <p 
in order to calculate the Hubble expansion rate as the 
function of the same parameters. We will then confront 
the type la supernova (SNIa), Baryon Acoustic Oscilla- 
tion (BAO) and Cosmic Microwave Background (CMB) 
data with the predictions of the model. 



II. UNPARTICLE SCALAR FIELD 

Consider a scalar unparticle degree of freedom 0, with 
the standard kinetic term and simple potential V(4>) ~ 
(f> 2 . The energy density and pressure can be written as 
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where M is the mass of this field. Since unparticles do 
not have a fixed mass, we take the mass M to be the 
average of the mass distribution 



J A " p{p)P 2 dp 2 
Jo" P(p) d P 2 



M 



(3) 



where p{p) is the mass spectral density of unparticles, 
while A u is the scale of the phase transition above which 
unparticle description stops being valid and therefore is 
the maximum mass of unparticles. It is important to 
distinguish X u from Mp. While Mp is the new funda- 
mental scale with ordinary particle degree of freedom, 
A„ is the scale at which these fundamental degrees of 
freedom are seen as unparticles. The normalization fac- 
tor is introduced to preserve the unparticle number. We 
believe that taking the average of the (continuous) un- 
particle mass distribution is the most convenient way to 
capture an average collective influence of unparticles on 
the evolution of the energy density of the universe. Of 
course, this is not a unique way to do this, one may for 
example consider two fields, one at the top and one at 
the bottom of the mass distribution. It will be interest- 
ing to consider alternative approaches, however, we take 
here the minimalistic average behavior approach. 

The mass spectral density of unparticles with continu- 
ous mass 1491 can be written as 



p(p) 



2</„ 



A 2d u 



(4) 



where d u is the scaling dimension of unparticles. The 
scale A„ is included for dimensional reasons and in the 



simplest case is equal to the "phase transition scale" X u . 
Strictly speaking the exact normalization of unparticles 
is unknown since we do not know how many unparticle 
degrees of freedom there are, neither we know the precise 
combinatorics factors (see discussion in [III). However, 
since we are taking an average of the mass distribution 
of unparticles in Eq. ([3]), precise normalization of unpar- 
ticles does not matter. Substituting (0| in ^ we get 
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We consider now an unparticle scalar field with negli- 
gible kinetic energy, which is an analog of a slowly rolling 
scalar field. From Eq. (p]), the energy density is 



Pu 



d u - 1 
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A 2 2 



(0) 



Now, we make the assumption that, as long as the 
kinetic energy is negligible, A„ w 0. This is just a state- 
ment that the value of the field is of the same order of 
magnitude as the maximum unparticle mass for the rele- 
vant energy scale. In this case we can write equation 
as 
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We can now express the mass term as 

1 



\m 2 = %— 

2 2d u 



(7) 



(8) 



The constant B is, in general, a function of the un- 
particle scaling dimension d u . If we set the constant of 
proportionality between A M and <j> to unity, we obtain 



B 



d u -l 
2d u 



(9) 



While this may not be an exact expression for B, it 
should be correct up to an 0(1) factor. 

The energy density and pressure can now be replaced 

by 



Pu 



dt 



By/p~u~(, 



Pu = -B^ 2 



(10) 
(11) 



In what follows, we will use this unparticle setup in the 
cosmological context. We will assume that the scalar un- 
particle degree of freedom plays the role of a quintessence 
field. 

We note here that in [5(| a different equation of state 
was derived for unparticle degrees of freedom. However, 
there, a model of bosons in thermal equilibrium was used, 
which is appropriate for the dark matter description, but 
not for our purpose of a slowly rolling scalar field. 
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III. THE MODEL 



We first rewrite the relation Eq. ([10]) for our purpose 
as (dots denote time derivatives): 



1 



Pv 



PA 



(12) 



where we have shifted the potential by the cosmological 
constant p\ for convenience. From here, one can express 
Pu as 



Pu 



B 2 (f> 4 + 20 2 + 4p A 



(13) 



The potential function of the scalar field V 
therefore 



is 



V 



BjB 2 fa + 2(/> 2 +4p A )(), 



The Lagrangian of the unparticle field is: 
1 



C 



( W - v 



T PA 

(14) 



(15) 



From the Euler-Lagrange equations, we obtain the fol- 
lowing equation of motion for <fi (where a subscripted 
variable denotes a partial derivative with respect to that 
variable): 



l-Vi 



3H - U 



V tt> -3HV i = (16) 



Now consider a Universe consisting of perfect fluid dark 
matter, radiation and unparticle dark energy. The dark 
matter energy density p m and the radiation density p r 
satisfy the continuity equations 



3Hp. n 







AHp r = 



(17) 



(18) 



The evolution of the expansion rate H = - (where a is 
the scale factor) is given by the Friedmann equation 



— (pm + ltf + V^) 



(19) 



The system of equations (|16p -(|T9 j) is solved numeri- 
cally. Initial conditions are set at some arbitrary rcdshift 
zt deep within the radiation dominated era. For the nu- 
merical analysis we work in units where c = K = 8irG = 
1, and we scale our time coordinate such that the present 
day Hubble value Hq = 1 for a AC DM cosmology (in 
what follows, 0-subscripts will always denote present day 
values). The initial value of the matter and radiation 
densities can be obtained by appropriately rcdshifting 
the observed present-day matter and radiation densities 
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FIG. 1: 1 + w(a) for different choices of B. fa is fixed to 1 



to Zj. The initial velocity of the field is set to zero, since 
the field can be assumed to be frozen on account of the 
very high Hubble friction - we have checked that even if 
the field is given an initial velocity it quickly comes to 
rest. The evolution is continued until the matter density 
parameter reaches the observed value of £l m o ~ 0.3. 

Clearly, the three parameters that completely deter- 
mine the evolution of the system are the following: 

• The observed present matter density parameter 

• The constant B, which is a function of the scaling 
dimension of unparticles, as explained above (see 
Eq. © and the discussion below it). 

• The initial value of the field fa. As per our assump- 
tions that <f> ~ A u , if the kinetic term is negligible, 
and that the field is slowly rolling, the parameter 
fa is essentially the same as initial X u . 

We now consider the evolution of the equation of state 
of the unparticle dark energy. The effect of varying B on 
the function w(a) is shown in Fig. |T]), where the initial 
value of the field is fixed at fa = 1. Note that w always 
starts out at —1, as a result of our choosing the initial 
velocity of the field to be zero. For small values of B, 
w increases slowly. However, as B increases, w evolves 
faster. If one chooses a smaller fixed value of fa, the plots 
have similar shapes, but smaller amplitudes. 

Fig. @ shows (on a logarithmic scale) the effect of 
varying the initial condition on (f> for a given choice of 
B = 0.25. Note that in our units, <j> is measured in units 
of the Planck Mass Mpi. We find that small values of fa 
lead to a very slow evolution of w, whereas large values 
lead to a stronger fluctuation in w. If one makes the same 
plots for smaller (or larger) fixed values of B, the shapes 
of the plots remain the same, though the amplitude of 
fluctuation of w is smaller (or larger). 

As explained in the next section, the parameter choices 
(B,fa) = (.5,1), correspond to the upper limits of the 
ranges of B and fa explored in this paper. From Fig. ([T]) 
we find that the final 1 + w in this extreme case is only 
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about 0.06. This shows that our slow-roll approximation 
is valid in all cases examined. 



IV. COMPARISON TO OBSERVATIONS 

In this section, we use observational data from Type la 
supernovae, Baryon Acoustic Oscillations and the Cosmic 
Microwave Background to constrain the three parameters 
il m0 , B and 6 t . 

a. Type la Supernovae constraints To place stan- 
dard candle constraints on the model, we use the recently 
released Union08 compilation of Snla data [l[ ■ This is a 
heterogeneous dataset consisting of data from the Su- 
pernova Legacy Survey, the Essence survey, the recently 
extended dataset of distant supernovae observed with the 
Hubble Space Telescope, as well as older datascts. 

The x 2 from SNIa is calculated as follows: 



N 5 
I] [Mobs (Zi) - /i t h (Zi)Y 

2 _ i=l 

XSN — 2 



(20) 



where N = 307 is the number of SNIa data points. /x OD s is 
the observed distance modulus, defined as the difference 
between the apparent and absolute magnitude of the su- 
pernova. The <7 M; i are the errors in the observed distance 
moduli, arising from a variety of sources, and assumed to 
be gaussian and uncorrclatcd. The theoretical distance 
modulus fi t h depends on the model parameters via the 
dimcnsionlcss luminosity distance D^(z): 

D l (z) = (l + z) f dz'——^—— (21) 



as follows: 
Mth (z) 



42.38 - 51og 10 ^ + 51og 10 [^ L (z)] (22) 



We marginalize over the present value of the Hubble 
parameter following the techniques described in (5l| and 
construct x 2 likelihood contours for the parameters f2 m o 
and B for different choices of Si. 



b. CMB constraints We use the CMB data to place 
constraints on the parameter space following the recipe 
described in [H]. The "CMB shift parameters" [H, [54[ 
arc defined as follows: 



R = (0)H o r , l a = nr( Z!f ) /r s (z r ) (23) 

Here r(z) is the comoving distance to redshift z defined 
as: 



r{z) 



H(z) 



dz 



(24) 



r s (z») is the comoving sound horizon at decoupling (red- 
shift z») given by 



0*) 



1 



H(z) ^3(1 + ^/(1 + 2)) 



zdz 



(25) 



The quantity Rj, is the photon-baryon energy-density 
ratio, and its value can be calculated as Rb = 
31500fl b h 2 (T CM B/2.7Ky 4 . The redshift at decoupling 
z* (fib, h) can be calculated from the formulas in [55j |. 

R can be physically interpreted as a scaled distance to 
recombination, and l a can be interpreted as the angular 
scale of the sound horizon at recombination. 

The x 2 contribution of the CMB is given by 



Xcmb 



Vcmb C inv VcMB 



(26) 



Here V C mb = P-Pdata, where P is the vector (l a ,R, z„) 
and the vector Pdata is formed from the WMAP 5-year 
maximum likelihood values of these quantities [12]. The 
inverse covariance matrix Ci nv is also provided in [52j ■ 

c. Baryon Acoustic Oscillation constraints Finally, 
we also use the Baryon Acoustic Oscillations (BAO) to 
constrain the model. The measured quantity here is the 
ratio r s (z*) jDy (z), where Dy (z) is the so called "vol- 
ume distance" defined in terms of the angular diameter 
distance Da = r (z) / (1 + z) as 



D v (z) 



(l + z) 2 D\(z)z 



-I 1/3 



H(z) 



(27) 



So far the BAO peak has been measured at two rcdshifts, 
z = 0.2 and z = 0.35 The ratio of the two measure- 

ments of D v (z), i.e., D v (.35) /D v (.2) = 1.812+0.060 [|, 
can be used as a model-independent observational con- 
straint. In this paper, following [5l|, we calculate the x 2 
contribution of the BAO measurements as follows: 



XbAO — V BAO C inv VbAO 



(28) 



The vector Vbao = P 
(D v (0.32) , D v (0.2)) and 
two measured BAO data points 
ance matrix is provided in Q. 



Pdata, with P = 
(0.1980,0.1094), the 
. The inverse covari- 
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FIG. 3: Likelihood plot from SNIa+BAO data for the pa- 
rameters B and J1a with cj>i = 0.01. The yellow (light) region 
is excluded at the 2<r level, and the darker (orange) region is 
excluded at the la level. Red (darkest) region is not excluded 
at either confidence level. 
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FIG. 4: As for Fig. © with fa = 0.5. 



FIG. 5: As for Fig. © with fa = 1 
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FIG. 6: Likelihood plot from CMB data for the parameters B 
and J1a with fa = 0.01. The yellow (light) region is excluded 
at the 2cr level, and the darker (orange) region is excluded at 
the lcr level. Red (darkest) region is not excluded at either 
confidence level. 



d. Constraint contours In figures [SHE] we construct 
joint SNIa+BAO constraint contours for the variables 
f^o and B for three different choices of fa, marginalizing 
over the present day Hubble parameter h. 

For the parameter B, we examine the range < B < 
0.5. This is in accordance with Eq. (|9|), using the fact 
that 1 < d u < 2 [28|. However, it is important to re- 
member that given our assumptions, Eq. ([9]) is true only 
up to an overall 0(1) factor. The limiting case of B = 
corresponds to a ACDM cosmology. 

In figures we construct (n,po,B) contours for three 
different values of fa (or equivalently, A), fa = .01,0.5, 1 
using the CMB data. Recall that in our units, <fi is mea- 
sured in units of Mpi, and hence the extreme case of 
fa = 1 corresponds to the maximum unparticlc mass be- 
ing on order of the Planck mass. 

Clearly, the observational data do not rule out these 
models. We find that the contours slightly shrink as 
fa increases. However, B is poorly constrained by the 
Snla+BAO data for all choices of fa examined, as well 
as the CMB data for small fa. This is because a smaller 
choice of fa leads to a slower evolution of the field for 



all choices of B (within the range examined). For large 
values of fa (closer to the Planck Mass), the CMB data 
place stringent constraints on B. 
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FIG. 7: As for Fig. (® with fa = 0.5 
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FIG. 8: As for Fig. © with fa = 1 



whose mass is protected. [We note however that a proper 
analysis of quantum corrections in the unparticle con- 
text would be highly useful.] Further, coupling of un- 
particles to the standard model fields can be arbitrar- 
ily suppressed by raising the fundamental energy scale 
Mp. Therefore, the unparticle dark energy model is free 
of most of the problems that plague fundamental scalar 
field quintessence models. 

We examined how the dynamics of the unparticle field 
equation of state is determined by the parameters Sl mj 0j 
the constant B (a function of the unparticle scaling di- 
mension d u ), and the initial value of the field 4n (which 
we assume to be comparable to the unparticle energy 
scale X u ). We have used the SNIa, CMB and BAO to 
place observational constraints on these parameters and 
find that such models are not ruled out by the data. 



CONCLUSIONS 



We have studied model of dark energy in which an un- 
particle degree of freedom plays the role of a slowly rolling 
scalar field quintessence. While fundamental scalar fields 
are not protected from radiative corrections, unparticles 
are, since they are bound states of fundamental fcrmions 
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